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Abstract 

Parity violating effects for low energy elastic neutron deuteron scattering are calculated for 
DDH and EFT-type of weak potentials in a Distorted Wave Born Approximation, using realistic 
hadronic strong interaction wave functions, obtained by solving three-body Faddeev equations in 
configuration space. The results of relation between physical observables and low energy constants 
can be used to fix low energy constants from experiments. Potential model dependencies of parity 
violating effects are discussed. 
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I. INTRODUCTION 



The study of parity violating (PV) effects in low energy physics is a very sensitive tool 
to test methods of calculations both of weak and strong interactions in the Standard model. 
This also can be a way to search for a possible manifestations of new physics resulted in 
deviations from unambiguous and precise calculations of PV effects and experimental mea- 
surements. However, to use this approach, it is crucial to prove that implemented theoretical 
techniques are sufficient to describe experimental data with high accuracy which exceeds ex- 
perimental accuracy. There is a large amount of experimental data for different PV effects in 
nuclear physics, each of which in general agrees with theoretical predictions. However, in the 
last years it became clear (see, for example l|-|4| and references therein) that the traditional 
DDH [5| method for calculation of PV effects cannot reliably describe the whole available set 
of experimental data within the same set of parameters. If this is not the manifestation of 
new physics, which is very unlikely for the current accuracy of experimental measurements 
and theoretical calculations, then this discrepancy could be blamed on systematic errors 
in experimental data, theoretical uncertainties in calculations of strong interactions at low 
energy, or it might be that DDH approach is not adequate for the description of the set of 
precise experimental data because it is based on a number of models and assumptions. To 
resolve this discrepancy and to eliminate nuclear model dependent factors in calculations, it 
is necessary to focus on the analysis of new and existing experimental data for different PV 
parameters in few-body systems, where calculations of nuclear related effects can be done 
with a high precision. Recently new approach, based on the effective field theory (EFT), 
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Deen introduced for a model independent parametrization of PV effects (see, papers 
4| and references therein), and some calculations for two-body systems have been done 
6| . The power of the EFT approach for parametrization of all P V effects in terms of a small 
number of constants could be utilized if we can analyze a large enough number of PV effects 
to be able to constrain all free parameters of the theory which are usually called low energy 
constants (LEC). Thus, one can guarantee the adequate description (parametrization) of 
the strong interaction hadronic parts and weak interaction constants for symmetry violat- 
ing observables. Unfortunately, the number of experimentally measured (and independent 
in terms of unknown LECs) PV effects in two body systems is not enough to constrain all 
LECs. In spite of the fact that five independent observable parameters in a two body system 
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could fix five unknown PV LECs 



10] , it is impossible to measure all of them using exist- 



ing experimental techniques. Therefore, one has to include into analysis few-body systems 
and even heavier nuclei, the latter of which are actually preferable from the experimental 



point of view, because as a rule, the measured effects in nuc 
nucleon-nucleon system due to nuclear enhancement factors 11 



ei are much larger than in 
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The natural and unambiguous way to verify the applicability of the EFT for the calcula- 
tion of symmetry violating effects in nuclear reactions requires a development of a regular 
and self consistent approach for calculation of PV amplitudes in three-body (few-body) sys- 



tems 



14(, with a hope to extend the formalism for the description of many body systems. 



This systematic approach for the solution of three-body PV scattering problem in EFT 
framework [l^] requires additional numerical efforts and will be presented elsewhere. As a 
first step for the clarification of the possible difference in contributions to PV effects from 
DDH and EFT-type potentials, one can use a "hybrid" method (similar to the method used 
in paper [15[) for the simplest process of neutron-deuteron scattering. We calculate three- 
body wave functions with realistic Hamiltonians of strong interaction using exact Faddeev 
equations in configuration space, and then, calculate PV effects in the first order of pertur- 
bation with DDH potential and potentials derived in EFT formalism. In the next section, 
we present our formalism for the calculation PV effects for elastic neutron-deuteron scatter- 
ing with different set of nucleon weak potentials, with DDH and weak potentials obtained 
from pionless and pionful EFTs. Then, we present results of numerical calculations and 
discussions. 



II. FORMALISM 

We treat weak nucleon interactions as a perturbation and calculate three-body wave func- 
tions exactly using Faddeev equations with phenomenological potentials for strong interac- 
tions. Similar hybrid approach has been successfully applied to the weak and electromagnetic 



processes involving three-body and four-body hadronic systems [16H21|. We consider three 
types of parity violating potentials. The first one is the standard DDH potential which 
is based on meson exchange mechanism of nucleon-nucleon interactions. The second and 
third potentials are derived from pionless and pionful versions of effective field theory with 
parity violating hadronic interactions. Instead of calculating parity violating amplitudes by 
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summing PV diagrams in EFT, we use these potentials to calculate PV effects. This is a 
simplification, which we call a "hybrid" approach. 



A. Observables 

Since PV effects in neutron-deuteron system are very small, we consider only coherent 
processes which are related to the propagation of neutrons through unpolarized deuteron 
target and, therefore, do not have an additional suppression in low energy region. Then, two 
PV observable parameters are the angle of rotation of neutron polarization around neutron 
momentum and the relative difference of total cross sections P = (er+ — er_)/ (er + + er_) for 
neutrons with opposite helicities. The value of the angle of neutron spin rotation per unit 
length of the target sample can be expressed in terms of elastic scattering amplitudes at 
zero angle for opposite helicities /+ and /_ as 

d<f) 2irN 

Tz =~— Re(/ + -/_), (1) 

where N is a number of target nuclei per unit volume and p is a relative neutron momentum. 
Using optical theorem, one can write the relative difference of total cross sections P in terms 
of these amplitudes as 

= Im (/+-/-) 

M/+ + /-)' {) 

It is convenient to represent the amplitudes in terms of matrix R which is related to 
scattering matrix S as R = 1 — S. With partial waves decomposition for the case of 
neutron-deuteron scattering 

\p,m n ,m d ) = J2 E \pAlyS)JJl(JJ%i;,SM)(SM\^m n Am d )Y* lZy (p), (3) 

lyl* SM,JJ* 

where l y is an orbital angular momentum between neutron and deuteron, S is a sum of 
neutron spin and deuteron total angular momentum, and J is the total angular momentum 
of the neutron-deuteron system, the above equations can be written at low energies as 

wi = Km + 4,j - - ^ii 
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and 
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u 2' 2 2' u 2 u 2' 2J 



/Re Ql + 2i2j a a , (5) 

u 2' u 2 u 2' u 2- 



where Ri> S i i S = (l'S'\R J \lS) , unprimed and primed parameters correspond to initial and 
final states. Since we are interested in low energy neutron scattering, it would be sufficient 
to include only s- and p-waves contributions to parity violating amplitudes; for the total 
cross section (the denominator in the last equation), we keep only dominant contributions 
from s-wave neutrons. It should be noted that time-reversal invariance leads to the relation 
(1S'\R J \0S) = (0S'\R J \1S') between matrix elements, therefore, only half of parity violating 
amplitudes are independent. 

Nucleon-nucleon interaction can be written as a sum V = V pc + V pv of the parity con- 
serving (Vp C ) and weak parity violating (V pv ) terms. Due to the weakness of parity violating 
interaction, one can use Distorted Wave Born Approximation (DWBA) to calculate PV 
amplitudes with a high level of accuracy as 

Ri>s>,i v s * 4r^+Vp h<*, (i' yS ')jr\v pv \y, {i y s)jr)£\ (6) 

where \i is a neutron-deuteron reduced mass and (l' y S')JJ z ) P c are solutions of 3-body 
Faddeev equations in configuration space for parity conserving strong interaction Hamilto- 
nian, defined by Vpc and normalized as described in section Hi CI The factor i~ l 'v +l y in this 
expression is introduced to match the i?-matrix definition in the modified spherical harmon- 
ics convention 22] with the wave functions which are calculated in this paper using spherical 
harmonics convention. 

In the rest of the paper, we use only wave functions calculated for parity conserving 
potentials and, therefore, will omit subscript PC. 

As will be explained in section Hi C\ we use jj-coupling scheme (with a basis states \l y j y )) 
when solving Faddeev equtions. One can transform jj-basis states into Z^iS-basis by means 
of 



\[ly® (S* ®jx)s]jj t ) = Yl ^ x ® ( l V ® S ^3y]jJz) 
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-l)i-+i«-J(_l)'v+«fc+j-+J[(2j 1/ + l)(2S + 1)]5 { ly Sk h \ ,(7) 

jx J S 



One interesting observation is that the neutron spin rotation, as well as parameter P, in 
\lyj y ) basis involves potential matrix elements only between j y = | states. 

It should be noted that at low energy the Im(i# 5 ,, s ) ~ p y + y+ , and thus the expres- 
sion eq.fll]) for the angle of neutron spin rotation is finite and well defined in the zero 
energy limit of the n-d scattering. Numerically, it is calculated by evaluating expression 
lm(Rf, s , lyS ) /p l 'y +ly+1 at zero energy. On the other hand, Re(i?^ 5 , ~ p ■ lm(Rf, s , lyS ) at 
low energy, and thus the real part of this quantity vanishes in the zero energy limit. There- 
fore, the parameter P is calculated at 15 KeV neutron kinetic energy in the laboratory 
system, where both imaginary and real parts of the R-matrix elements become comparable 
in magnitude and thus can be discerned numerically. 



B. The parity violating potentials 

To understand the possible difference in the description of parity violating effects by DDH 
and EFT-type for potentials, we compare calculations with the DDH potential|5] and two 



different choices of EFT potentials: the potential derived from pionless E 



T lagrangian [l| 

and the potential derived from pionful EFT Lagrangian It was shown isj] that all these 
three potentials can be expanded in terms of a set of Off operators as 



>.i 



c n°ij\ a = DDH or pionless EFT or pionful EFT (8) 



with parameters c" and operators given in the Table HI 

One can see that operators O^j are products of isospin, spin, and vector operators X^ ± 
defined as 

= [Py,/n(^)] + , 



-(») 



Xll'- = ilPijJnirij)}-, (9) 



where p {j — 

For the DDH potential, radial functions f x {r), x = n,p, and u are modified Yukawa 
functions, 



A r r ( mi 



(10) 
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TABLE I: Parameters and operators of parity violating potentials. irNN coupling g^NN can be 
represented by qa by using Goldberger-Treiman relation, g n = gAfn^/F^ with F n = 92.4 MeV. 
Tij = (3t?t? — n ■ Tj). Scalar function L\(r) = 3La(v) — H^{r). 
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For pionless EFT (/EFT) one, /„(r) are described by single function / M (r), 

1 



Ur) 



Airr 



(11) 



with /i ~ m^. 

For the case of pionful EFT model (7rEFT), there are long range interactions from one 
pion exchange(y_i ; £#) and from their corrections (Vi : lr), middle range interactions due to 
two pion exchange (V^mr), an d short range interactions (V^sr) due to nucleon contact 
terms. The radial part of the leading term of long range one pion exchange, V-^lr , is 
described by the function Ur). Since one-pion exchange contribution is dominated by long 
range part, we do not use a regulator for it, i.e. we assume that the long range interactions 
have the same radial functions Ur) as DDH potential with infinite cutoff. The short range 



interaction V\ t sR in pionful theory has the same structure as for pionless EFT; however, 
in spite of the structural similarity, their meanings are rather different. One can ignore 
the higher order corrections of long range interactions, V^lr, because they can either be 
absorbed by renormalization of low energy constants [6] or suppressed. The middle range 
interactions V^mr are described by functions L(q) and H(q) in momentum space 

^.v^JvMlZiM), g(g)s ±f L(g)> (12) 

W \q\ \ 2m, y W Aml + q 2 W ' V ; 

where, = (q°, q) = Px—Pi = V2 ~ P%- To calculate two pion exchange functions (divergent 
at large q) in spacial representation , we use regulators Tpq^pp ■ For the sake of simplicity, 
we use only one cutoff parameter with the same regulator, both for middle range and for 
short range interactions. Then, one can write 

{L A (r), H A (r), / A (r)} = ± J ^ e -^ ^l^£ {L( g ), H{q), 1}. (13) 

In the given representation, coefficients c" have fm dimension and scalar functions f%(r) have 
fm _1 dimension. One can see that only the new operator structure, which is not included 
in DDH or pionless EFT, is due to V\lr- Therefore, pionful EFT does not introduce new 
operator structure, provided we neglect V^r term [6|? ]. 

To see a sensitivity to the choice of cutoffs for parity violating potentials, we used two 
set of cutoff parameters for each models, which are listed in the Table HI1 

Using the discussed above three potentials, one can represent parity violating amplitudes 
as a linear expansion in terms of given set of matrix elements for corresponding operators 
O^. Thus, the angle of neutron spin rotation can be written as 

i* = f cT (14) 



and the parameter P as 



Ndz 

n=l 



13 



p = ^«, ( 15 ) 



n=l 



in terms of coefficients J£ and i£ with a = DDH-I,II, /rEFT-1,11, ttEFT-1,11 for different 
potentials and cutoff parameters. 
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TABLE II: Parameter of parity violating potentials in GeV units. We used masses of mesons m. 
m p , and m u , respectively, as 0.138, 0.771, and 0.783 in DDH potential. 







K 






n 




A 


DDH-I 


1.72 


1.31 


1.50 


>EFT-I 


0.138 


ttEFT-1 


0.8 


DDH-II 


oo 


oo 


oo 


>EFT-II 


1.0 


ttEFT-II 


1.0 



C. Faddeev wave function 



To obtain 3-body wave functions for neutron-deuteron scattering with parity conserving 



interactions, we solve 
configuration space 



23 



addeev equations (also often called Kowalski-Noyes equations) in 



241 ] . For isospin invariant interactions (with nucleon masses fixed 
to H 2 /m = 41.471 MeV-fm), three Faddeev equations become formally identical, having the 
form 

(E-Hq- Vn) = MA + (16) 
where (ijk) are particle indices, Hq is kinetic energy operator, is two body force between 
particles i, and j, ipk = "4>ij,k is Faddeev component. 

The wave function in Faddeev formalism is the sum of three Faddeev components, 

*(aj, y) = ipxixx, y x ) + ^2(^2, y 2 ) + ^3(^3, y 3 )- (17) 

Using relative Jacobi coordinates x k = (rj — r^) and y k = -^(^k ~ r '^ rj )? one can expand 
these Faddeev components in bipolar harmonic basis: 

rl>k = E Fa{Xk,Vk) (k (s iSj ) s ) . (l y s k ) ) ® t k ) TT , (18) 

where index a represents all allowed combinations of the quantum numbers presented in 
the brackets: l x and l y are the partial angular momenta associated with respective Jacobi 
coordinates, Sj and t% are the spins and isospins of the individual particles. Functions 
F a (x k ,y k ) are called partial Faddeev amplitudes. It should be noted that the total angular 
momentum J as well as its projection M are conserved, but the total isospin T of the system 
is not conserved due to the presence of charge dependent terms in nuclear interactions. 

Boundary conditions for Eq. (fT6|) can be written in the Dirichlet form. Thus, Faddeev 
amplitudes satisfy the regularity conditions: 

F a (0,y k ) = F a {x k ,0) = 0. (19) 



'ilk 



For neutron-deuteron scattering with energies below the break-up threshold, Faddeev com- 
ponents vanish for x& — > oo. If yfc — > oo, then interactions between the particle k and the 
cluster ij are negligible, and Faddeev components ipi an d 4>j vanish. Then, for the compo- 
nent ipk, which describes the plane wave of the particle k with respect to the bound particle 
pair ij, 

lim ^ fc (x fc ,y fc ) /nin = —= Y]|{0d(x fc )}. ® {^(y*) ® ) ® {Utj) t t^) 

X ^ [^-Wn^W ~ SlUMnK^nd) > (2°) 

where deuteron, being formed from nucleons z and j, has quantum numbers = 1, jd — 1, 
and = 0, and its wave function 0d(xfc) is normalized to unity. Here, r n d = (v3/2)yfc 
is relative distance between neutron and deuteron target, and h t are the spherical Hankel 
functions. The expression ( 120]) is normalized to satisfy a condition of unit flux for nd 
scattering wave function. 

For the cases where Urbana type three-nucleon interaction (TNI) is included, we modify 
the Faddeev equation (1TB]) into 

(E-H - V v ) ^ = Vytyi + V»i) + \(Vj k + (21) 



by noting that the TNI among particles ijk can be written as sum of three terms: Vi 



v)k 



D. Evaluation of matrix elements 



Due to anti-symmetry of the total wave function in isospin basis, one has (\I/|Vi2 + V23 + 
V31I*) = 3(^1^1^) for any pair i^j. 
Using decomposition of momentum p, 



p = -i\/ x = -i x— + -Vn 

OX X 



(22) 



we can represent general matrix elements of local two-body parity violating potential oper- 
ators as 



V3, 



<-><*/ioi*«> (+) = (Y v ' 



a/3 



dxx dyy 



X x) 



xy 



xy 



(a\0(x)\m3) 
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where (±) means outgoing and incoming boundary conditions and X(x) is derivative of 
scalar function or derivative of wave function with respect to x. (Note that we have used 
the fact that (F^)* = F^ + \) The partial amplitudes Fj(/) jQ ,(x, y) represent the total systems 
wave function in one selected basis set among three possible angular momentum coupling 
sequences for three particle angular momenta: 

The "angular" part of the matrix element is 

(a\6(x)\P) = fdx f dyyt(x,y)6(x)yp(x,y), (25) 



where y a {x,y) is a tensor bipolar spherical harmonic with a quantum number a. One can 
see that operators for "angular" matrix elements have the following structure: 

O(x) = (r 2 r 3 )(cr 2 © o- 3 ) • x , or (r 2 t 3 )(<t 2 © cr 3 ) ■ V n , (26) 

where 0,® = ±, x. The explicit values of these matrix elements are summarized in the 
appendix. 

III. RESULTS AND DISCUSSIONS 

As it was mentioned in the previous section, because of low energy property of R J a i a1 it 
is convenient to present results for elements Rf, ly in terms of a ratio, 

" — H <*, (Zp)JJ M \VT\% (l y S)JJ z )^ (27) 



For the case of parity violation, we fix V — 1 and l y = 0. To obtain the observable 
parameters when neutron energies are larger than thermal ones (which correspond to zero 
energy limit for neutron spin rotation), one can use a simple extrapolation based on the 
above representation with a good accuracy up to hundreds KeV. 



The contributions to parity violating matrix elements -— Im 



R J , (p) 

a' a Ki 

2~ 



! from different terms 
of parity violating potentials (see Table [I]) are presented in the Table III II These matrix ele- 
ments were calculated using strong AV18+UIX and weak DDH-II parity violating potentials 
for the case of low neutron energies (up to thermal ones). From this table, one can see that 
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rJ is', jM 



at very low energy in fm 2 units. We chose AV18+UIX 



TABLE III: Contributions to -Im 

7T 

as strong potential and DDH-II as parity violating potential. Matrix elements of n = 6, 7 are zero 
due to of isospin structure. 



n 
1 

2 
3 
4 
5 
8 
9 

10 
11 
12 
13 



<?' - 1 j - I 

u — 2 , o — 2 



<?' - 3 T - I 
o — 2 ; <J — 2 



<?' - 1 f - 3 

o — 2 ' ^ — 2 



c/ _ 3 t _ 3 
o — 2 , «y — 2 



-0.522 X 10+ 00 

0.480 x 10~ 03 
-0.284 x 10~ 03 
0.797 x 10~ 04 
0.664 x 10~ 04 
-0.228 x 10~ 03 
0.278 x 10~ 04 
0.755 x 10~ 04 
0.631 x 10~ 04 
-0.105 x 10~ 01 
-0.175 x 10~ 01 



0.253 x 10+ 00 
-0.182 x 10~ 01 
0.339 x 10~ 02 
0.410 x 10~ 02 
0.475 x 10~ 02 
0.190 x 10~ 02 
-0.562 x 10~ 02 
0.388 x 10~ 02 
0.453 x 10~ 02 
0.452 x 10~ 02 
0.725 x 10~ 02 



0.131 x 10+ 00 
-0.105 x 10~ 01 
0.231 x 10~ 01 
-0.154 x 10~ 01 
-0.178 x 10~ 01 
0.180 x 10~ 01 
0.960 x 10~ 02 
-0.146 x 10~ 01 
-0.170 x 10~ 01 
0.165 x 10" 03 
0.113 x 10~ 02 



-0.151 x 10~ 01 
0.882 x 10~ 02 

-0.428 x 10~ 02 
0.221 x 10" 03 
0.313 x 10~ 03 

-0.301 x 10~ 02 
0.107 x 10~ 02 
0.209 x 10~ 03 
0.298 x 10~ 03 

-0.223 x 10~ 03 

-0.377 x 10-° 3 



the main contribution to PV effects comes from J = 3/2 channel for the "best values" of 
DDH coupling constants. 

Our results for the angle of neutron spin rotation for DDH, pionless EFT, and pionful 
EFT weak interaction potentials with different sets of parameters are summarized in Tables 
ITvWl an d lYU For these calculations, we used two types of strong interacting potentials: 
Argonne two nucleon interaction AV18 and inclusion of Urbana IX three nucleon interaction, 
AV18+UIX. One can see that these results practically do not depend on a choice of the strong 
interaction potential. Also, it is clear that the matrix element related to pion-exchange 
(n = 1) is dominant for DDH potential, slightly enhanced for pionfull potential, and about 
equal to other ones for pionless potential. 

The neutron spin asymmetry P was calculated for laboratory neutron energy E = 15 
KeV. The results are summarized in tables IVHl IVIIIt and IIXI for DDH, pionless EFT, and 
pionful EFT weak interaction potentials with different sets of parameters, correspondingly. 
These results provide a pattern similar to that of the results for the angle of neutron spin 
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TABLE IV: Coefficients 1% DH for AV18 and AV18+UIX strong potentials, and DDH-I and DDH-II 
parameter sets for parity violating potentials. = 0. 



n 


DDH-I/AV18 


DDH-I/ AV18+UIX 


DDH-II/AV18 DDH-II/AV18+UIX 


1 


0.612 x 10+ 02 




0.596 x 10+ 02 


0.616 x 10+ 02 


0.600 x 10+ 02 


2 


0.666 x 10+ 00 




0.726 x 10+ 00 


0.114 x 10+ 01 


0.124 x 10+ 01 


3 


-0.130 x 10+ 01 




0.133 x 10+ 01 


-0.212 x 10+ 01 


-0.217 x 10+ 01 


4 


0.911 x 10+ 00 




0.934 x 10+ 00 


0.131 x 10+ 01 


0.134 x 10+ 01 


5 


0.980 x 10+ 00 




0.992 x 10+ 00 


0.153 x 10+ 01 


0.156 x 10+ 01 


8 


-0.125 x 10+ 01 




0.130 x 10+ 01 


-0.160 x 10+ 01 


-0.167 x 10+ 01 


9 


-0.615 x 10+ 00 




0.622 x 10+ 00 


-0.786 x 10+ 00 


-0.796 x 10+ 00 


10 


0.998 x 10+ 00 




0.102 x 10+ 01 


0.124 x 10+ 01 


0.127 x 10+ 01 


11 


0.111 x 10+ 01 




0.113 x 10+ 01 


0.146 x 10+ 01 


0.149 x 10+ 01 


12 


0.991 x 10+ 00 




0.983 x 10+ 00 


0.141 x 10+ 01 


0.140 x 10+ 01 


13 


0.144 x 10+ 01 




0.144 x 10+ 01 


0.226 x 10+ 01 


0.225 x 10+ 01 


TABLE V: Coefficients H 


for AV18 and AV18+UIX strong potentials, and >EFT-I and /rEFT-II 


parameter sets for parity violating potentials. 5 6 7 10 n 12 13 = 0- 




n 


/rEFT-I/AV18 


>EFT-I/AV18+UIX 


AEFT-II/AV18 /rEFT-II/AV18+UIX 


1 


0.616 x 10+ 02 


0.600 x 10+ 02 


0.969 x 10+ 00 


0.969 x 10+ 00 


4 


0.606 x 10+ 02 


0.588 x 10+ 02 


0.499 x 10+ 00 


0.515 x 10+ 00 


8 


-0.761 x 10+ 02 


-0.757 x 10+ 02 


-0.677 x 10+ 00 


-0.708 x 10+ 00 


9 


-0.946 x 10+ 01 


-0.662 x 10+ 01 


-0.341 x 10+ 00 


-0.348 x 10+ 00 


rotation. The parameter J n in these tables is defined as 






J n = Re 

C n TV 


4/xp 2 


I ol - 2V2RI o1 

l 2 > 2 2 ' 2 


+ 44 o3 -2y/ERh Q3 ) 

1 2' 2 2 ' 2 / 


(28) 


and 


is related to the parameter I 


n in the expression P = c nln by 










§(27r/V)J n 


_ 87T 2 /i J n 


(29) 






/„ = 

Re 


- 1. 3 

u 2' 2 2 ' 2 . 


9 (J tot ' 



13 



TABLE VI: Coefficients I£ for AV18 and AV18+UIX strong potentials, and vrEFT-I and vrEFT-II 
parameter sets for parity violating potentials. 3 6 7 10 11 12 = 0- 



n 


TTtjb 1-1/AVlo 


vrEFT-I/AV18+UIX 


7TEFT-II/AV18 


^_T7 T7T^ TT / A^\7"1 O i T TT V 

TTShr 1-11/AVlo+UlA 


1 


0.616 x 10 +u 


0.600 x 10 H 


-02 


0.616 x 10 H 


-02 


0.600 x 10 +u 


A 

4 


0.152 x 10 + 


n 1 a \y 1 rH 
U.14z X 1U 


-01 


A KA O \y 1 rH 

U.04y X 1U 


-00 


r\ A ci -I 1— III) 

0.488 x 10 + 


5 


0.435 x 10+ 01 


0.185 x 10 H 


-01 


0.123 x 10 H 


-01 


0.664 x 10 -01 


8 


-0.184 x 10+ 01 


-0.179 x 10 H 


-01 


-0.782 x 10 H 


-00 


-0.748 x 10+ 00 


9 


-0.820 x 10+ 00 


-0.730 x 10 H 


-00 


-0.340 x 10 H 


-00 


-0.288 x 10+ 00 


13 


0.226 x 10+ 02 


0.218 x 10 H 


-02 


0.970 x 10 H 


-01 


0.936 x 10+ 01 


14 


0.339 x 10+ 01 


0.333 x 10 H 


-01 


0.177 x 10 H 


-01 


0.174 x 10+ 01 


15 


0.654 x 10+ 02 


0.631 x 10 H 


-02 


0.273 x 10 H 


-02 


0.264 x 10+ 02 



where Otot is the total n — d cross section. The total cross section o t ot can be calculated, or 
one can use its known experimental value. 

From the presented data, one can see that the results of our calculations are only slightly 
different for the cases when we use AV18 and AV18+UIX strong Hamiltonians. This indi- 
cates stability of the results with respect to the three nucleon forces. Indeed, by analyzing 
the DDH one-pion exchange matrix element (see Table IIIII) . one can see that for DDH-I 
with potentials AV18 and AV18+UIX, the contributions to the I n= \ are —0.180 x 10 +01 and 
—0.333 x 10 +01 for doublet channel (J = 1/2), and for the quartet channel (J = 3/2) they 
are 0.630 x 10 +02 and 0.630 x 10 +02 , correspondingly. The quartet channel is dominated by 
the repulsive and long-range part of the strong interactions, but the doublet channel is de- 
fined by attractive part. Therefore the quartet channel is less sensitive to the off-energy shell 
structure of the strong interactions compared to the doublet channel. Then, due to the dom- 
inant contribution from the quartet channel, the net result turns to be rather independent 
on the contribution from three nucleon forces. This fact demonstrates the independence of 
our results on models of strong interactions. However, further investigations with different 
strong interaction potentials are desirable. 

It should be noted, that the dependence on cutoff parameters for the contributions from 
potentials with short and middle range interactions, even though it appears large, does not 
lead to cutoff dependence for the observable parameters. Indeed, the renormalization of 
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TABLE VII: Coefficients J» DH for AV18 and AV18+UIX strong potentials, and DDH-I and DDH- 
II parameter sets for parity violating potentials at E = 15 KeV in the laboratory frame. J§ ® H = 0. 



n DDH-I/AV18 DDH-I/ AV18+UIX DDH-II/AV18 DDH-II/AV18+UIX 



1 


0.253 


X 


10 +oo 


0.253 


X 


10+ 00 


0.254 


X 


10+ 00 


0.254 


X 


10+ 00 


2 


0.246 


X 


1Q -02 


0.245 


X 


10 -02 


0.390 


X 


1Q -02 


0.384 


X 


1Q -02 


3 


-0.190 


X 


10 -02 


-0.147 


X 


10 -02 


-0.313 


X 


1Q -02 


-0.243 


X 


1Q -02 


4 


0.769 


X 


1Q -03 


0.393 


X 


10 -03 


0.110 


X 


10 -02 


0.563 


X 


10 -03 


5 


0.846 


X 


1Q -03 


0.442 


X 


1Q -03 


0.132 


X 


1Q -02 


0.689 


X 


1Q -03 


8 


-0.176 


X 


1Q -02 


-0.134 


X 


1Q -02 


-0.228 


X 


1Q -02 


-0.175 


X 


1Q -02 


9 


-0.235 


X 


1Q -03 


0.567 


X 


1Q -04 


-0.259 


X 


1Q -03 


0.118 


X 


1Q -03 


10 


0.842 


X 


1Q -03 


0.430 


X 


1Q -03 


0.104 


X 


1Q -02 


0.534 


X 


1Q -03 


11 


0.957 


X 


1Q -03 


0.500 


X 


1Q -03 


0.126 


X 


1Q -02 


0.657 


X 


1Q -03 


12 


0.374 


X 


1Q -02 


0.370 


X 


1Q -02 


0.528 


X 


1Q -02 


0.522 


X 


1Q -02 


13 


0.563 


X 


1Q -02 


0.559 


X 


1Q -02 


0.874 


X 


1Q -02 


0.868 


X 


1Q -02 



TABLE VIII: Coefficients jf t for AV18 and AV18+UIX strong potentials, and >EFT-I and >EFT- 
II parameter sets for parity violating potentials. Jf^ 5 6 7 10 n 12 13 = 



n /rEFT-I/AV18 

1 0.254 x 10+ 00 

4 0.503 x 10~ 01 

8 -0.111 x 10+ 00 

9 -0.241 x 10~ 02 



/rEFT-I/AV18+UIX 
0.254 x 10+ 00 
0.240 x 10~ 01 
-0.854 x 10- 01 
0.338 x 10~ 02 



>EFT-II/AV18 
0.372 x 10-° 2 
0.421 x 10-° 3 
-0.984 x 10~ 03 
-0.904 x 10~ 04 



>EFT-II/AV18+UIX 
0.369 x 10~ 02 
0.215 x 10-° 3 
-0.763 x 10~ 03 
0.750 x 10~ 04 



low energy constants would cancel those cutoff dependencies by the cutoff dependencies of 
LECs. Therefore, as a result, calculated PV observables are practically cutoff independent. 

All these tables present information about contributions of different PV operators to PV 
effects, provided we know corresponding weak coupling constants. Then, to calculate parity 
violating effects, we can use either DDH potential or one of the considered EFT potentials. 
However, for the case of EFT potentials, we need to know a set of LECs which cannot be 
calculated in the given theoretical framework but must be obtained from a number of inde- 
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TABLE IX: Coefficients J% for AV18 and AV18+UIX strong potentials, and vrEFT-I and ttEFT-II 
parameter sets for parity violating potentials. JJ 3 6 7 10 11 12 = 0- 



n 


ttEFT-I/AV18 


vrEFT-I/AV18+UIX 


7TEFT-II/AV18 


7TEFT-II/AV18+UIX 


1 


0.254 x 10+ 00 


0.254 x 10+ 00 


0.254 x 10+ 00 


0.254 x 10+ 00 


4 


0.106 x 10~ 02 


0.352 x 10~ 03 


0.309 x 10~ 03 


0.333 x 10~ 04 


5 


0.741 x 10" 02 


0.512 x 10~ 02 


0.292 x 10~ 02 


0.221 x 10~ 02 


8 


-0.276 x 10~ 02 


-0.212 x 10~ 02 


-0.127 x 10~ 02 


-0.100 x 10~ 02 


9 


-0.148 x 10~ 03 


0.301 x 10" 03 


-0.278 x 10~ 04 


0.168 x 10~ 03 


13 


0.976 x 10~ 01 


0.981 x 10~ 01 


0.421 x 10~ 01 


0.423 x 10~ 01 


14 


0.137 x 10~ 01 


0.136 x 10~ 01 


0.714 x 10~ 02 


0.712 x 10~ 02 


15 


0.283 x 10+ 00 


0.284 x 10+ 00 


0.119 x 10+ 00 


0.120 x 10+ 00 



pendent experiments. Unfortunately, currently available experimental data are not enough 
;o define the LECs with required precision. Even for pionless EFT, the estimated LECs 
1] have large uncertainties preventing us from predicting the values of PV effects. For the 
pionful EFT, the situation with determination of LECs is even worse. Therefore, it is impos- 
sible to make reliable predictions for PV effects using EFT-type potentials at this time, and 
the only reasonable way to estimate magnitudes of PV effects is to use the DDH potential. 
Taking into account the difficulty of the systematic description of PV effects using "stan- 
dard" DDH potentials (see discussions in the introduction), we estimate PV effects using 
the DDH potential for different sets of weak coupling constants: both for the "best value" 
coupling constants and for two possible sets of the values of the coupling constants recently 



obtained by Bowman 25] from the fit of reliable existing experimental data (see Table IXl) . 
The results for these three sets of weak coupling constants are summarized in Tables IXll and 
IXHI for the angle of spin rotation and for neutron spin asymmetry, correspondingly. One can 
see that in contrast to the fact that the one-pion exchange dominates in the DDH- "best" 
coupling parameter set, the rho meson exchange dominates in the case of Bowman's coupling 
parameter set. One can see that the angle of neutron spin rotation has almost the same 
magnitude for all three sets of parameters, but it has opposite signs for the "best value" 
set and for the Bowman's fits. The neutron spin asymmetry does not only have opposite 
signs but also essentially different values for these two choices of parameters. This allows 
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TABLE X: DDH PV coupling constants in units of 10 7 . Strong couplings are % = 13.9, % = 0.84, 

2 

!^ = 20, k p = 3.7, and k u = 0, h' p contribution is neglected. 4-paramter fir and 3-parameter fit 
uses the same h l p and hi with DDH 'best'. 

DDH Coupling DDH 'best' 4-parameter fit [25] 3-parameter fit[25] 

hi +4.56 -0.456 -0.5 

h° p -11.4 -43.3 -33 

h 2 p -9.5 37.1 41 

hi -1.9 13.7 

hp -0.19 -0.19 -0.19 

hh -1.14 -1.14 -1.14 



TABLE XI: Neutron spin rotation in 10 rad-cm for the case of DDH-II potential with 
AV18+UIX strong potential for a liquid deuteron density N = 0.4 x 10 23 atoms per cm 3 . 





DDH 


'best' 


4-parameter 


fit[25j 


3-parameter 


fit [25] 


1 


0.108 


X 


10+ 00 


-0 


108 


X 


10 -oi 


-0 


118 


X 


10 -oi 


2 


0.386 


X 


10 -02 





147 


X 


io- 01 





112 


X 


io- 01 


3 


-0.317 


X 


io- 01 


-0 


120 


X 


10 +oo 


-0 


918 


X 


io- 01 


4 


0.349 


X 


1Q -04 





349 


X 


1Q -04 





349 


X 


1Q -04 


5 


0.150 


X 


10 -03 





150 


X 


1Q -03 





150 


X 


1Q -03 


8 


-0.423 


X 


10 -02 





305 


X 


io- 01 





000 


X 


10 +oo 


9 


-0.202 


X 


1Q -02 





146 


X 


io- 01 





000 


X 


10 +oo 


10 


0.967 


X 


1Q -03 





967 


X 


1Q -03 





967 


X 


1Q -03 


11 


0.113 


X 


10 -02 





113 


X 


10 -02 





113 


X 


1Q -02 


12 


0.102 


X 


10 -02 





102 


X 


1Q -02 





102 


X 


1Q -02 


total 


0.768 


X 


10 -oi 


-0 


682 


X 


10 -oi 


-0 


891 


X 


10 -oi 



one to choose between two possible sets of DDH parameters and, as a consequence, to test 
the dominance of pion-meson contribution in PV effects in n — d scattering. 

Finally, we would like to mention that our results are quite different from the results 
obtained in paper 15] . For example, in paper 15) , the values of I n for J = \ and J = | have 
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TABLE XII: Neutron spin asymmetry for the case of DDH-II potential with AV18+UIX strong 



potential (the 


total cross section a tot 


= o.oo b at hi = Lt> Jvev J. 






DDH 'best' 


4- parameter fit [25] 


3- parameter fit [251 


1 


0.947 x 10-° 8 


-0.947 x 10~ 09 


-0.104 x 10~ 08 


2 


0.248 x 10~ 09 


0.943 x 10~ 09 


0.719 x 10~ 09 


3 


-0.740 x 10~ 09 


-0.281 x 10~ 08 


-0.214 x 10~ 08 


4 


0.304 x 10~ 12 


0.304 x 10~ 12 


0.304 x 10~ 12 


5 


0.138 x lO" 11 


0.138 x lO" 11 


0.138 x lO -11 


8 


-0.922 x 10~ 10 


0.665 x 10~ 09 


0.000 x 10 +0 ° 


9 


0.620 x lO" 11 


-0.447 x 10~ 10 


-0.000 x 10+ 00 


10 


0.843 x lO" 11 


0.843 x 10~ n 


0.843 x 10~ n 


11 


0.104 x 10~ 10 


0.104 x 10~ 10 


0.104 x 10~ 10 


12 


0.797 x 10~ 10 


0.797 x 10~ 10 


0.797 x 10~ 10 


total 


0.899 x 10~ 08 


-0.209 x 10~ 08 


-0.236 x 10~ 08 



the same signs for operator with n — 1, but our results show opposite signs for these matrix 
elements. Another discrepancy is related to the systematic difference between the values of 
matrix elements calculated [15J for AV18 and AV18+UIX potentials, which indicates a large 
wave function difference for AV18 and AV18+UIX potentials. Contrary to those, our results 
show that these matrix elements are insensitive to the presence of the three nucleon force. 1 



IV. CONCLUSION 



We have calculated parity violating angle of neutron spin rotation and asymmetry in 
transmission of neutrons with opposite helicities for low energy neutron deuteron scattering. 
Using Distorted Wave Born Approximation for weak interactions with realistic three-nucleon 
wave functions from Faddeev equations in configuration space, we have parameterized PV 
observables in terms of matrix elements presented in the DDH weak potential and in weak 
potentials derived from pionless and pionful EFTs. It is shown that our results practically 



1 We thank R. Schiavilla and M. Viviani for discussions which clarified that the reason for theses discrep- 



ancies is related to numerical errors in paper 15]. 
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do not depend on the choice for strong interaction potentials and on cutoff parameters. 

Based on the given analysis, one can see that for DDH potential, the dominant contri- 
bution to observable PV effects comes from the pion-exchange matrix element with n — 1. 
However, for pionless EFT potential, all types of matrix elements contribute almost equally, 
and for pionful EFT potential the pion-exchange matrix element is sightly enhanced as 
compared to the other ones. Therefore, it would be interesting to compare the estimation 
of observable PV effects using appropriate LECs and coupling constants for DDH. Un- 
fortunately, due to insufficient data for LECs this is impossible at this time. However, a 
comparison of PV effects for two different sets of coupling constants shows that n — d scatter- 
ing experimental results can be used to distinguish between different sets of DDH coupling 
constants and to help in clarification of the issue about the importance of the contribution 
of pion-exchange weak potential. 

Appendix: Explicit results of angular part of Matrix elements 

Explicit values of matrix elements of iso-spin operators for two-body states are 



and matrix elements of orbital and spin operators for two-body states \(l x s x )jxjx) are 



(TX\n-T 2 \TT z ) 
(T%\( ri + T 2 y\TT z ) 

(t't' z \(ti-t 2 ) z \tt z ) 

(T'T'Mn x t 2 ) z \TT z ) 



^t z ,t z ^t',t[^t,i — 35t,o], 
8t>\8ti$t' z t z [2T z ], 
5t',t±iSt z ,t z St z ,o[^\, 
£ T , Tz S Tz ,oSt' ,t±i [±2] , 




(A.l) 



({jx ± 1, IjWi^i + 0-2) • x\(Jx, l)j x f x ) 



(Ux, l)jx£\(<Tl + <T2) • X\(jx ± 1, l)jxjx) 




(A.2) 



(Ux ± 1, - ^2) • X\(jx, 0)jx3x) 

((jx, 0)j x f x \((T 1 - (T 2 ) ■ X\(j x ± 1, l)j x f x ) 




(A.3) 
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((jx, 0)j x f x \i((T 1 X <T 2 ) • x\(j x ± 1, l)j x j x ) 

= (-)((jx ± 1, l)j x 3 z x \i(<Ti x <J 2 ) ■ x\(J x , 0)j x f x ) 

= ±2 J^ + l ' 2±l ' 2 (A.4) 

(om^i(^^ 2 )-v,jo, ± n M ) = , 2 ^« 

^ imv-n (j, + 1/2 T 1/2)^, + 1/2 ± 1/2 
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